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ABSTRACT 

We present new estimates of the spatial two-point correlation function of rich clusters of 
galaxies selected from the APM Galaxy Survey. We have measured redshifts for a sample 
of 364 clusters out to a depth of ~ 450 /}~^Mpc*. The clusters have a mean space density 
of n = 3.4 X 10~® /i^Mpc" . The two-point correlation function, ^cc, for this sample is 
equal to unity at a pair-separation of ro = 14.3 ± 1.75 /}~"'^Mpc (2(7 errors), consistent 
with our earlier results from a smaller sample. The new observations provide an accurate 
determination of the shape of ^^c to pair-separations of about 50 /}~^Mpc. Our results 
show that ^cc has a higher amplitude than expected according to the standard il = 1 
cold dark matter (CDM) model on spatial scales 2 < s < 50 /}~^Mpc, but are in 
good agreement with scale-invariant fluctuations in either a low density CDM model 
or a critical density universe made up of a mixture of hot and cold dark matter. Our 
results provide strong constraints on so called 'co-operative' models of galaxy formation 
in which the galaxy formation process introduces large-scale structure in the galaxy 
distribution. 
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INTRODUCTION 



Rich clusters of galaxies have been used as tracers of large- 
scale structure by a number of authors over the last few 
years. Many of these studies have been based on Abell's 
(1958) cluster catalogue, which was compiled by visual in- 
spection of uncalibrated photographic plates. However, there 
is considerable evidence that the clustering measured from 
Abell's catalogue is affected by inhomogeneities in the se- 
lection of clusters on the plane of the sky (Sutherland 1988, 
Soltan 1988, Efstathiou et al. 1992). As a result, several 
groups have attempted to produce more uniform samples 
by applying computer algorithms to select rich clusters of 
galaxies from calibrated photographic plates (Dalton et al. 
1992, hereafter DEMS, Nichol et al. 1992) or by construct- 
ing X-ray flux limited samples (Romer et al. 1993, Nichol, 
Briel & Henry 1994). The results from these newer surveys 
give 



i.e. the correlation length, ro, is lower than the value ro ^ 
20 h~ Mpc deduced from redshift surveys of Abell clusters 
(Bahcall & Soneira 1983, Klypin & Kopylov 1983, Postman, 
Huchra & Geller 1992), though Bahcall & West (1992) have 
argued that the discrepancy may be a consequence of a cor- 
relation between ro and cluster richness rather than inhomo- 
geneities in the Abell catalogue. The computer-selected and 
X-ray selected redshift surveys have so far been limited to 
small samples of typically 100-200 clusters and so the shape 
of S,cc is relatively poorly constrained, especially on scales 
> 15 A" Mpc where i^cc is smaller than unity. 

In this Short Communication we present the results of 
a new redshift survey of 364 rich clusters of galaxies selected 
from the APM Galaxy Survey (Maddox et al. 1990a, Mad- 
dox, Efstathiou & Sutherland 1990). This is a extension of 
our earlier redshift survey of 220 APM clusters described 
in DEMS and Dalton et al. (1994). The larger sample de- 
scribed here allows us to determine the shape of the S^cc more 
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Where h is Hubble's constant in units of lOOkni s 



Mpc" 
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accurately and to larger scales than in our earlier work. The 
space density of the cluster sample described here is compa- 
rable to that of our original survey (DEMS) and so in this 
paper we concentrate on a comparison of S^cc with theoreti- 
cal predictions rather than investigating the dependence of 
S,cc with cluster richness. The richness dependence of S^cc will 
be discussed in a future paper using a new redshifts for a 
deeper sample of particularly rich APM clusters. 

The new redshift survey is described in Section 2, and 
estimates of the two-point correlation function are presented 
in Section 3. In Section 4 we compare our results with the 
predictions of various CDM-like models, using N-body simu- 
lations as described by Croft & Efstathiou (1994a). We show 
that the new results allow us to draw strong conclusions on 
the parameters of acceptable CDM-like models. Our results 
also provide constraints on models where a physical process 
other than gravity modulates the efficiency of galaxy forma- 
tion on large scales, and so changes the large-scale properties 
of the galaxy distribution. (Rees 1986, Babul & White 1991, 
Bower et al. 1993). 



2 THE CLUSTER SAMPLES 

DEMS have described an algorithm for selecting rich clusters 
of galaxies from the APM Survey. Here we describe some 
changes to the DEMS algorithm that were made to produce 
an approximately volume limited cluster sample to greater 
depths so allowing us to extend the DEMS redshift survey. 

The depth of the cluster sample used by DEMS was 
limited by the magnitude range used to define the cluster 
richness. The cluster richness, 7?., was defined by DEMS to 
be the weighted number of galaxies in the magnitude range 
{mx — 0.5, mx + 1.0) above the mean background count in 
the range {mx — 0.5, rrex + 1.5), where mx is the magni- 
tude of the galaxy for which the weighted count above back- 
ground exceeded X = 7?./2. The cluster catalogue defined in 
this way is limited by the point at which the faint end of the 
background slice equals the magnitude limit of the survey, 
which occurs when m,x = 19.0. We increased the depth of 
our new cluster catalogue by changing the background slice 
to (mx — 0.5, mx + 1.0), and redefining X to be TZ/2.V As 
a starting point for our new catalogue we used the same set 
of percolation centres as DEMS. A detailed description of 
the cluster finding algorithm and a discussion of the effects 
of changing the various selection parameters will be given in 
a future paper (see also Dalton 1992). 

In the new cluster catalogue 117 clusters with 7?. > 50 
and mx < 19.2 had redshifts measured from our earlier red- 
shift survey, Dalton et al. (1994). Dalton et al. (1994) list 
redshifts for an additional 20 clusters which have 7?. > 50 
and 19.2 < mx < 19.5. This left 223 clusters without red- 
shifts to form a complete sample of clusters with 7?. > 50 and 
mx < 19.2. We observed typically three galaxies in the cen- 
tral regions of 154 clusters in four nights using the 4m tele- 
scope at Cerro Tololo Inter- American Observatory (CTIO). 



T We shall use the symbols mx and 7Z to refer to magnitude 
and richness estimates within the catalogue defined here; hence- 
forth we will use the subscript DEMS to refer to the analogous 
quantities defined for the DEMS catalogue. 
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Figure 1. The redshift distributions of the two cluster sam- 
ples discussed in the text. The solid line in each panel shows 
the smoothed distribution (see text). The dot-dashed lines show 
estimates of the selection functions (arbitrary normalisations) ob- 
tained by dividing the smoothed distribution by the volume ele- 
ment for an Einstein-de Sitter universe. 



We applied the likelihood ratio test of Dalton et al. (1994) 
to each redshift measurement. We were unable to determine 
the redshift of 26 clusters because of poor quality spectra, 
or discrepancies in the galaxy redshifts none of which sat- 
isfied the likelihood ratio test as probable cluster members. 
We shall include the 20 clusters with previously measured 
redshifts, but with mx > 19.2, and refer to this sample of 
265 clusters as sample A. The mean space-density measured 
for this sample using the method described in Section 2.3 of 
Efstathiou et al. (1992) is 2.2 x 10"^ /^^Mpc"^ similar to 
the space density of 2.4 x 10~ h Mpc~ estimated for the 
DEMS TIdems > 20 sample by Efstathiou et al. (1992). 

We can increase the size of the redshift sample by trans- 
forming the richnesses of clusters in the DEMS sample with 
measured redshifts into the richness TZ defined here. Com- 
paring the richnesses of clusters measured in both cata- 
logues, we find the relation 



7^ = 24.8 -I- 1.27^I 



(2) 



with a scatter of 16 in TZ. Transforming all of the DEMS 
richnesses to TZ using equation 2 and selecting clusters with 
7?. > 50 leads to a sample of 364 clusters with measured red- 
shifts, which we call sample B. The mean space density of 
clusters in sample B is approximately 3.4 x 10~ h Mpc~ . 
This sample is not limited strictly according to TZ as defined 
by the cluster selection. However it defines a statistical sam- 
ple that is almost equivalent to a 7?. > 50 sample, because the 
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Figure 2. The correlation functions for samples A and B. The 
filled circles show the estimates obtained using Hamilton's (1993) 
estimator. The power law model ^cc = {r/14 /i~ Mpc)~ which 
provided a good fit to the results of our earlier survey (Dalton 
et al. 1992) is shown for comparison. 
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Figure 3. The correlation function Ccc{o': ^) for sample B shown 
as a function of pair separations perpendicular to the line of sight 
((t) and along the line of sight (tt). Contour levels are at ^cc = 
3, 2, 1,0. 8, 0.6, 0.4, 0.2,0. ,-0.2, -0.4. The contour level at (cc = 
1 is shown by the heavy line; negative contours are plotted as 
dashed lines. 



scatter in the transformed richnesses is similar to the scatter 
intrinsic to the cluster finding algorithm (Dalton 1992). 



3 THE CLUSTER CORRELATION FUNCTION 

We estimate the redshift-space correlation functions for the 
two samples using the estimator 



Uis) = 2/ 
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(DR) 



1, 



(3) 



where / is the ratio of the number of random points to the 
number of clusters in the sample. In each case we use 20, 000 
points distributed within the survey boundaries and with 
the same redshift distributions as the smoothed distribu- 
tions shown in Figure 1 which were obtained by convolving 
the observed redshift histograms with Gaussians of width 
4000km s~ . We estimated ^cc using clusters with redshifts 
in the range in the range 5000-50000km s~ . The correlation 
functions for the two samples are shown in Figure 2, together 
with estimates obtained from sample B using Hamilton's ( 
1993) estimator: 
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(DR)'' 

which is less affected by uncertainties in the selection func- 
tion for S^cc < 1. In equations 3 and 4 the numerical factor 
accounts for the fact that we only count each (DD) or (RR) 
pair once. The error bars plotted in th e Figur e are computed 
from the formula S^cc = (1 + icc)/\/(DD) which Croft & 
Efstathiou (1994a) have shown underestimates the variances 
in simulations of rich cluster catalogues in CDM-like models 
by a factor of between 1.3 and 1.7. 
Fitting a power law 



Uis) 



to the estimates of ^cc over the range 1 h ^Mpc ^ s < 
40 h~ Mpc for samples A and B plotted in Figure 2 gives: 



7 = 2.14 ±0.19, 
7 = 2.05 ±0.12, 
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where the error bars give an approximate indication of the 
95% confidence intervals on the parameters of the fits. The 
agreement between samples A and B and the excellent agree- 
ment of the fits with those quoted by Efstathiou et al. (1992) 
shows that the correlation function is insensitive to changes 
in the cluster finding algorithm. 

Figure 3 shows contour plots of the the correlation func- 
tion for Sample B measured as a function of the pair sep- 
aration along the line of sight, tt, and perpendicular to the 
line of sight, a. This Figure shows that S,cc{o', tt) is nearly 
isotropic. There is no evidence of the large anisotropics mea- 
sured in redshift surveys of Abell clusters by Efstathiou et al. 
(1992) in which correlations of order unity are measured to 
distances of tt ~ 100 h~ Mpc. The estimates of ^cc{o', tt) 
for sample A are almost identical to those plotted in Fig- 
ure 3 for sample B. The absence of anisotropics in the APM 
cluster catalogue provides strong evidence that the machine 
selected cluster samples are free of inhomogeneities on the 
plane of the sky. 



4 COMPARISON WITH THEORY AND 
DISCUSSION 

In Figure 4 we compare the estimates of S^cc for sample 
B with theoretical predictions from dissipationless N-body 
simulations of CDM-like models (Croft & Efstathiou 1994a, 
Croft & Efstathiou 1994b). The three models shown are: (a) 
standard CDM, i.e. scale-invariant initial flucuations in an 
critical density universe with h = 0.5; (b) alow density CDM 
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Figure 4. The data for sample B (filled circles in figure 2) compared with the average correlation functions of clusters selected from 
ten N-body simulations each of (a) Standard CDM, (b) a low density CDM model with non-zero cosmological constant, and (c) a mixed 
dark matter model (see text). 



model (LCDM) with Q = 0.2, h = 1 and non-zero cosmo- 
logical constant introduced to make the universe spatially 
flat; (c) a mixed dark matter (MDM) model with h = 0.5 in 
which CDM contributes Qcdm = 0.6, light neutrinos con- 
tribute Q^ = 0.3 and baryons contribute Q_b = 0.1. The 
details of the simulations are given in Croft & Efstathiou 
(1994a) and Croft & Efstathiou (1994b). 

We locate rich clusters in real space in the simulations 
by first identifying cluster centres using a percolation algo- 
rithm and then computing the mass enclosed within a radius 
Tc = 0.5 h~ Mpc of the cluster centre-of-mass. The clus- 
ter catalogues are then ordered by mass and a lower mass 
limit is applied to the catalogue so that the mean space 
density of the clusters matches the observed space density 
(3.4 X 10"^ A^Mpc"^ for sample B). Croft & Efstathiou 
(1994a) describe the algorithm used to select clusters in 
more detail and show that the two-point cluster correlation 
functions measured from the simulations are insensitive to 
the cluster identification radius Tc, to the amplitude of the 
mass fluctuations and, at the space-densities of interest here, 
to cluster richness. 

Each panel in Figure 4 shows S^cc determined from the 
simulations for two values of the amplitude of mass fluc- 
tuations as measured by ag, which gives the rms of the 
mass density fluctuations 8p/p measured in spheres of ra- 
dius 8 h~^Mpc. As shown by Croft & Efstathiou (1994a), the 
shape of S^cc is remarkably insensitive to the value of ag . We 
note here that the COBE background radiation anisotropy 
measurements (Wright et al. 1994) require cg ~ 1 for mod- 
els (a) and (b) and ag Si 0.67 for model (c) (see Croft & 
Efstathiou 1994b) if the initial fluctuation spectrum is as- 
sumed to be scale invariant and temperature anisotropics 
from gravitational waves can be ignored. We have therefore 
used the simulation outputs for ag = 1 for models (a) and 
(b) and ag = 0.6 for model (c) to construct mock APM 
cluster catalogues in redshift space with the observed clus- 
ter redshift distribution and the APM Survey boundary. The 



two-point correlation functions measured from the simulated 
APM catalogues are plotted as the thick solid lines in Fig- 
ure 4. The redshift-space cluster correlation functions for 
models (b) and (c) are almost indistinguishable from the 
real-space estimates, except on scales > 30 h~ Mpc, where 
the redshift-space estimates are in better agreement with 
the observations. The difference between the redshift space 
and real space correlation function are more pronounced for 
model (a) with ag = 1, though even here they are still rela- 
tively small. 

Figure 4 shows that our new estimates of S^cc lie higher 
than the predictions of the standard CDM model by Si 3a on 
all scales plotted in the Figure (i.e. between 2 h~ Mpc and 
50 h-^Mpc). The LCDM and MDM models give a much 
better match to the observed cluster correlations. In fact, 
the MDM model provides a slightly better fit to the obser- 
vations than the Q/t = 0.2 LCDM results plotted in Figure 4. 
However, an LCDM model with a lower value of Q/t, may 
provide a more acceptable fit to the data (a numerical sim- 
ulation of LCDM with Q/t = 0.15 gave results in excellent 
agreement with the observations). 

The LCDM and MDM models provide a good match to 
the galaxy correlations on scales ^ 5 h~ Mpc measured from 
the APM Survey (Maddox et al. 1990b; Efstathiou, Bond & 
White 1992) if the galaxy two-point correlation function is 
related to the mass autocorrelation function by a linear bias- 



ing model, fg 
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vhere 6 is a constant. Some authors 



(e.g. Babul & White 1991) have argued that the galaxy for- 
mation process might introduce correlations on large scales, 
in which case the true relation between the galaxy and mass 
autocorrelation functions might be more complicated than 
assumed in the linear biasing model. For example, large- 
scale modulations in the galaxy distribution might arise if 
galaxies form with reduced or enhanced efficiency near to 
quasars (Rees 1986, Babul & White 1991, Efstathiou 1992, 
Bower et al. 1993). 

The results shown in Figure 4 provide strong evidence 
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against such models for the following reasons: (i) the cor- 
relation functions of rich clusters and galaxies have simi- 
lar shapes but different amplitudes, as expected if clusters 
and galaxies are linearly biased tracers of the mass distribu- 
tion; (ii) predictions of the rich cluster correlation function 
should be free of many of the uncertainties associated with 
the physics of galaxy formation if clusters form in regions 
of strong mass concentrations; (iii) the large differences be- 
tween the observed form of S^cc and the predictions of the 
standard CDM model, even on scales where S^cc is greater 
than unity would seem difficult to explain by a spatial mod- 
ulation of the efficiency of galaxy formation. 

In summary, we have extended the APM cluster redshift 
survey to 364 rich clusters, sampling an effective volume of 
Si fO h~ Mpc . The correlation function of rich clusters 
measured from this sample is in excellent agreement with 
our earlier results and is well approximated by a power law 
fee = {ro/sf with To = f4.3 ± f.75 h~^Mpc. N-body sim- 
ulations of CDM-like models show that the standard CDM 
model fails to match the observations by a wide margin. 
The data are in good agreement with an MDM model with 
^CDM ~ 0.6 and Q^ Si 0.3, or with a low-density variant of 
CDM with QA ~ 0.2. Croft & Efstathiou (f994b) have ar- 
gued that the peculiar velocities of clusters favour the MDM 
model over alow density CDM model, and analyses of galaxy 
streaming motions also favour Q = f universes (Kaiser et al. 
f99f, Dekel et al. f993). A consistent picture of structure 
formation is therefore beginning to emerge in which Q = f 
and the mass distribution is more strongly clustered than 
expected according to the standard CDM model. 
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